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1 . INTRODUCTION 


Nuclear  weapons  deployed  at  high  altitudes  produce  high-altitude 
electromagnetic  pulse  effects  (HEMP)  that  can  seriously  degrade 
tactical  weapon  and  communication  systems  vitally  needed  by  the  field 
army  prepared  to  fight  a conventional  and  nuclear  war.  The  Multiple 
Systems  Evaluation  Program  (MSEP)  was  established  to  determine  both  the 
vulnerability  and  the  means  for  hardening  a large  number  of  U.S.  Army 
tactical  systems  to  an  HEMP  environment.  An  essential  step  in  the 
program  is  to  develop  analytic  tools  to  evaluate  system  susceptibility 
to  HEMP,  that  is,  computer  programs  for  predicting  transient  data  and 
system  response.  These  computer  programs  have  been  gathered  into  an 
applications  package  titled  GAMPHS1  (Generic  Assessment  Methods  for 
a Priori  Hardening  of  Systems) . The  GAMPHS  application  for  the 
vulnerability  and  hardness  assessment  of  systems  covered  by  MSEP  uses 
two  modular  computer  libraries,  TRANL  and  ANAPAC,  in  addition  to  other 
computer  programs  (fig.  1) . These  library  routines  are  the  part  of 
GAMPHS,  whose  purpose  is  to  provide  proven  analytical  computer  programs 
to  commodity  commands  or  project  managers  for  assessment  of  the  EMP 
vulnerability  of  Army  systems.  The  application  of  TRANL  and  ANAPAC  to 
GAMPHS  is  explained  in  section  4 of - this  report. 

The  TRANL  and  ANAPAC  modular  analysis  package  libraries  are  a 
collection  of  independent  computer  software  programs  and  subroutines, 
each  written  to  do  a specific  task,  and  stored  in  two  permanent  files 
that  form  user's  libraries.  The  modular  building  blocks  in  ANAPAC  can 
be  combined  in  the  user's  computer  program  to  meet  his  needs.  Most  of 
the  computational  building  blocks  are  input/output-free  to  allow  the 
user  greater  flexibility  in  writing  a driver  program.  Two  of  the 
modular  subroutines,  RDTAPE  and  WRTAPE , are  provided  for  input/output 
compatibility  of  data  files  written  on  either  disc  or  tape.  Except  for 
the  fast  Fourier  transform  algorithm,  the  software  is  written  in 
FORTRAN  IV.  The  fast  Fourier  transform  algorithm  is  written  in  COMPASS 
assembly  language,  which  noticeably  decreases  its  execution  time. 

\ 

The  TRANL  library  is  a collection  of  independent  programs  that 
primarily  supports  the  data  digitizing  system.  The  programs  can  be 
merged  into  the  control  streams  of  major  programs  to  preprocess 
digitized  data  and  prepare  the  data  for  analysis.  The  data  structure 
of  the  TRANL  programs  is  compatible  with  RDTAPE  and  WRTAPE.  The  use  of 
programs  TEDD  and  REDUCE  is  not  restricted  to  digitized  data.  All  the 
programs  in  TRANL  are  written  in  FORTRAN  IV. 


1 George  Gornak  et  al,  EMP  Assessment  for  Army  Tactical  Communi- 
cations Systems:  Transmission  Systems  Series  No.  1 — Radio  Terminal  Set 

AN/TRC-145  (U) , Harry  Diamond  Laboratories  TR-1746  (February  1976). 
(SECRET -RESTRICTED  DATA) 
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MULTIPLE  SYSTEMS  EVALUATION  PROGRAM 

COMPUTER  PROGRAM  FLOW  IN  VULNERABILITY  ASSESSMENT 


Figure  1.  The  GAMPH3  application  in  vulnerability  and  hardness. 


The  ANAPAC  and  TRANL  libraries  are  stored  on  the  Control  Data 
Corporation  (CDC)  Computer  System  at  the  U.S.  Army  Mobility  Equipment 
Research  and  Development  Command  (USAMERADCOM)  and  accessed  by  means  of 
the  ATTACH  and  LIBRARY  Scope  control  cards.  The  subroutines  in  ANAPAC 
are  accessed  by  a standard  FORTRAN  CALL  statement.  The  programs  in 
TRANL  can  be  accessed  only  through  the  control  cards. 


2.  DESCRIPTION  OF  PROGRAMS  AVAILABLE  IN  THE  ANAPAC  LIBRARY 
2.1  Plotting  Routines 

The  plotting  software  presented  here  was  written  for  use  with 
a Houston  Instrument  DP-1  plotter  on  line  to  a Mohawk  2400  system.  The 
Mohawk  system  is  used  as  a remote  job  entry  (RJE)  station  to  a 
CDC  6600. 
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2.1.1  Subroutine  DRAW1 


The  FORTRAN  callable  computer  subprogram  DRAW12  is  designed 
to  plot  a single  graph  of  one  function  (data  set)  complete  with  axes, 
necessary  scaling,  and  user-defined  annotations  (axis  labels,  plot 
title,  and  plot  subtitle) . The  plot  can  be  a linear-linear  plot, 
log-log  plot,  linear-log  plot,  or  log-linear  plot.  When  arrays  are  to 
be  plotted  on  a log  axis,  the  software  will  drop  values  that  are 
invalid  for  a logarithmic  scale  (that  is,  zero  or  negative  values)  when 
they  occur  at  the  beginning  of  an  array;  a message  is  then  printed  that 
notifies  the  user  what  has  been  done.  Invalid  values  occurring  within 
an  array  cause  the  termination  of  the  plot,  and  an  appropriate  message 
is  printed.  The  date  and  time  of  execution  for  the  processing  of  each 
graph  are  displayed  in  the  lower  right-hand  corner  of  the  graph.  The 
first  time  DRAW1  is  called,  a remote  operator  message  (in  this  system, 
a "PM  message")  is  generated  that  alerts  the  operator  that  THIS  JOB 
PLOTS. 

User-defined  annotation  refers  to  axis  labels,  plot  title, 
and  plot  subtitle,  all  of  which  are  optional.  Axis  labels  may  contain 
up  to  three  words  (30  characters) . When  an  array  contains  axis  labels 
that  are  greater  than  three  words,  the  output  to  the  plotter  is 
unpredictable . Besides  the  user-defined  axis  notations,  the  plotting 

software  adds  a scaling  factor  after  any  label,  if  necessary,  to 
prevent  the  numerical  notations  on  the  axis  from  overlapping  illegibly. 
This  scaling  factor  notation  is  recorded  in  the  following  format: 

+N 

"(X  10  )" 

The  plot  title  and  subtitle  may  each  consist  of  up  to  eight 
words  (80  characters) , placed  at  the  top  of  the  plot.  The  subtitle  is 
placed  below  the  title  in  slightly  smaller  print.  To  center  the  title 
or  subtitle,  the  user  fills  the  appropriate  arrays  with  leading  blanks. 
The  user  can  incorporate  special  characters  (see  fig.  2)  into  the  title 
and  subtitle  by  punching  the  corresponding  Hollerith  punch  in  the 
appropriate  character  position  in  the  title  and  subtitle  arrays. 

For  the  standard  and  default  axis  position,  the  plotter 
places  the  two  axes  so  that  they  intersect  in  the  lower  left  corner  of 
the  graph.  For  arrays  having  both  positive  and  negative  values,  two 
other  options  are  available  for  the  linear  axis.  The  axis  can  be  drawn 
either  at  the  zero  point  occurrence  or  at  the  standard  location  with  a 
parallel  line  at  the  zero  occurrence  as  a reference  line.  If  these 
options  are  selected  and  the  array  values  do  not  span  zero,  the 
standard  axis  position  is  selected  by  default. 

2Thomas  V.  Noon,  Enhanced  Plotting  Software  for  Use  with  the  Houston 
Instrument  Complot  Plotter,  Harry  Diamond  Laboratories  TM-75-32 
(December  1975) . 
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Figure  2.  Hollerith  code  for  special  characters. 


If  the  number  of  points  to  be  plotted  is  passed  to  DRAW1  as 
a positive  value,  the  data  set  arrays  are  restored  to  the  values  they 
had  before  being  passed  to  DRAW1.  This  is  done  by  recalculating  the 
values.  If  the  values  of  the  data  set  arrays  are  no  longer  needed,  the 
number  of  points  to  be  plotted  can  be  passed  as  a negative  value.  This 
will  result  in  the  data  arrays  not  being  restored,  thus  saving  some 
execution  time. 
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The  first  six  formal  parameters  in  the  subroutine  argument 
list  are  integer  variables;  the  remaining  eight  are  real  variables. 
The  definition  and  use  of  each  is  as  follows: 

ITYPE  Variable  used  to  select  the  type  of  plot. 

If  ITYPE=1  a linear-linear  plot  is  selected. 

=2  a log-log  plot  is  selected. 

=3  a linear-log  plot  (semilog  with  the  linear 
axis  for  the  dependent  array)  is  selected. 

=4  a log-linear  plot  (semilog  with  the  linear 
axis  for  the  independent  array)  is  selected. 

IXLAB  Number  of  words  in  the  label  for  the  independent  axis,  not 
to  exceed  3.  If  IXLAB=0,  axis  label  is  omitted. 

IYLAB  Number  of  words  in  the  label  for  the  dependent  axis,  not  to 
exceed  3.  If  IYLAB=0,  axis  label  is  omitted. 

IPTLAB  Number  of  words  in  the  plot  title,  not  to  exceed  8.  If 
IPTLAB=0,  plot  title  is  omitted. 

ISBLAB  Number  of  words  in  the  plot  subtitle,  not  to  exceed  8.  If 
ISBLAB=0,  plot  subtitle  is  omitted. 

NPT  Number  of  points  to  be  plotted  in  both  the  independent  and 
dependent  data  arrays.  If  the  negative  of  the  number  of 
points  is  passed  through  NPT,  the  data  arrays  are  not 
restored  to  the  values  they  had  before  being  passed, 

ZAX  Variable  that  determines  where  the  independent  axis  is 
drawn . 

If  ZAX=0 . axis  is  drawn  in  the  standard  location.  This  is 
also  the  default  position  when  the  next  two 
conditions  cannot  be  met. 
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-1.  axis  is  drawn  at  the  zero  occurrence  of  the 
dependent  axis. 

=2,  axis  is  drawn  at  a standard  location  with  a zero 
reference  line  placed  at  the  zero  occurrence  on 
the  dependent  axis. 

Variable  that  determines  where  the  dependent  axis  is  drawn. 

If  ZAY=0.  axis  is  drawn  in  standard  location.  This  is 
also  the  default  position  when  the  next 
conditions  cannot  be  met. 


=1.  axis  is  drawn  at  the  zero  occurrence  of  inde- 
pendent axis. 

=2.  axis  is  drawn  at  standard  location  with  the 
addition  of  a zero  reference  line  placed  at  the 
zero  occurrence  on  the  independent  axis. 

Array  containing  user-defined  label  for  the  independent 
axis,  not  to  exceed  3 words  (30  characters) . If  IXLAB=0,  a 
dummy  name  can  be  used  in  the  call. 

Array  containing  user-defined  label  for  the  dependent  axis, 
not  to  exceed  3 words  (30  characters) . If  IYLAB=0,  a dummy 
name  can  be  used  in  the  call. 

Array  containing  user-defined  title  for  the  plot,  not  to 
exceed  8 words  (80  characters).  If  IPTLAB=0,  a dummy  name 
can  be  used  in  the  call. 

Array  containing  user-defined  subtitle  for  the  plot,  not  to 
exceed  8 words  (80  characters).  If  ISBLAB=0,  a dummy  name 
can  be  used  in  the  call. 


Array  containing  values  of  independent  variable  to  be 
plotted. 


The  following  informative  messages  are  printed  to  the  line 
printer  by  DRAW1: 


(a)  "XX  POINTS  DROPPED  FROM  BEGINNING  OF  TRACE,  NEGATIVE 
OR  ZERO  VALUES  IN  ARRAY" 


This  message  tells  the  user  that  there  is  some  number  (XX) 
of  invalid  values  (negative  or  zero)  at  the  beginning  of  an  array  for 
the  plotting  of  the  array  on  a log  scale.  The  graph  will  start  at  the 
first  valid  value. 

(b)  "TRACE  CANNOT  BE  PLOTTED,  NEGATIVE  OR  ZERO  VALUES  IN 
ARRAY" 

This  message  is  printed  when  an  invalid  value  (negative  or 
zero)  occurred  after  valid  values  for  the  plotting  of  the  array  on  a 
log  scale.  The  graph  will  not  be  plotted. 


(c)  "AMIN=XX  AMAX=YY" 


This  message  provides  the  user  with  additional  information 
on  the  minimum  and  maximum  values  of  an  array  when  their  difference  is 
very  small  and  expansion  of  the  resolution  of  the  scale  would  cause  the 
scale  notation  to  overlap  illegibly. 

2.1.2  Subroutine  DRAW4 


The  FORTRAN  callable  computer  subprogram  DRAW42  is 
designed  to  plot  one  or  more  functions  (data  sets)  on  the  same  graph. 
Any  number  of  data  sets  can  be  plotted  on  each  graph,  and  up  to  eight 
graphs  can  be  processed  "simultaneously."  The  axes,  scaling, 
user-defined  annotations  (axis  labels,  plot  title,  and  plot  subtitle), 
and  the  type  of  plots  available  (linear-linear,  log-log,  linear-log, 
and  log-linear)  are  processed  and  displayed  by  the  plotter  in  the  same 
manner  as  in  the  subprogram  DRAW1  (see  sect.  2.1.1).  As  in  DRAW1 , the 
date  and  time  of  execution  of  each  graph  are  displayed  in  the  lower 
right-hand  corner  of  the  graph,  and  the  first  time  the  subprogram  is 
called,  a PM  message  is  generated  to  indicate  that  "THIS  JOB  PLOTS." 
However,  DRAW4  has  greater  capabilities  than  DRAW1  for  the  presentation 
of  the  data  sets  and  the  appearance  of  the  overall  plot. 

2 Thomas  V.  Noon,  Enhanced  Plotting  Software  for  Use  with  the  Houston 
Instrument  Complot  Plotter,  Harry  Diamond  Laboratories  TM-7S-32 
(December  1975) . 
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To  provide  the  multiple  plotting  capability,  DRAW4  is 
divided  into  three  major  sections.  The  first  section,  accessed  by 
setting  the  first  formal  parameter  (IP)  to  1,  sets  up  the  user-defined 
annotations  and  a scratch  disc  file  identification.  This  section  also 
rewinds  and  initializes  the  scratch  file,  allowing  the  file  to  be 
reused  without  user  intervention.  The  second  section,  accessed  by 
setting  IP  to  2,  processes  each  data  set  and  the  manner  the  data  are  to 
be  plotted,  and  stores  this  information  on  the  scratch  disc  file 
initialized  in  the  first  section.  The  second  section  must  be  called 
for  each  data  set  to  be  plotted  and  can  be  called  any  number  of  times. 
The  third  section,  accessed  by  setting  IP  to  3,  completes  the 
processing  of  all  the  information  provided  in  the  first  and  second 
sections  for  the  specified  scratch  disc  file  and  sends  the  information 
to  the  OUTPUT  file. 

This  procedure  requires  N+2  calls  to  the  subprogram  DRAW4 , 
where  N is  the  number  of  data  sets  to  be  plotted.  Both  the  first  and 
third  sections  are  called  once  for  each  plot.  By  using  this  procedure, 
DRAW4  can  process  up  to  eight  unique  graphs  simultaneously,  thereby 
avoiding  the  need  to  recalculate  information  that  is  to  be  plotted 
later  along  with  other  information  in  the  program.  After  closing  a 
graph  by  a call  to  DRAW4  with  IP=3,  the  scratch  disc  file  used  for  that 
particular  graph  is  available  for  processing  another  graph  or  any  other 
application  required  in  the  program. 

The  first  section  of  DRAW4  (IP=1) , as  well  as  establishing 
the  user-defined  annotations,  declares  the  scratch  disc  file  to  be  used 
for  the  graph.  The  user  supplies  the  name  and  size  of  the  arrays  that 
contain  the  labels  of  the  independent  and  dependent  axes  and  the  plot 
title  and  subtitle.  Ii  a value  of  zero  is  passed  for  the  size  of  any 
one  of  the  above-mentioned  arrays,  that  notation  is  omitted  from  the 
graph.  The  label  information  supplied  is  then  stored  on  the  declared 
scratch  disc  file  for  use  in  the  third  section. 

The  second  section  of  DRAW4  (IP=2)  receives  the  data  set  to 
be  plotted  and  provides  the  user  with  a number  of  choices  for  the 
manner  in  which  a particular  data  set  is  displayed  on  the  graph.  A 
data  set  can  be  plotted  with  a solid  line,  a dashed  line,  a dotted 
line,  or  a dash-dot  line.  In  addition  to  these  line  types,  all 
standard  characters  and  various  symbols  (seven  are  provided)  can  be 
used  to  mark  the  occurrence  of  actual  data  points  in  a data  set.  These 
symbols  can  be  used  by  themselves  or  with  one  of  the  line  types.  The 
frequency  of  the  plotting  of  the  markers  can  also  be  varied  (e.g., 
every  data  point,  every  tenth  data  point,  etc).  The  second  section  of 
DRAW4  also  allows  for  the  specification  of  minimum  and  maximum  values 
for  the  dependent  array,  so  that  various  graphs  may  be  compared.  If 
the  specified  values  fall  outside  the  range  of  the  dependent  values  of 
all  the  data  sets  to  be  plotted,  these  specified  values  are  used  in  the 
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third  section  of  DRAW4  to  establish  the  scale  of  the  dependent  axis. 
If  the  specified  values  are  inside  the  range  of  the  dependent  arrays, 
these  arrays  are  changed  to  conform  with  the  given  minimum  or  maximum. 
Otherwise,  the  data  arrays  passed  to  this  section  are  not  altered. 

The  third  section  of  DRAW4  (IP=3)  allows  the  user  more 
control  over  the  drawing  of  the  axes,  zero  reference  lines,  and  border, 
and  over  the  selection  of  scale  factors  and  the  color  of  ink  used  to 
plot  a particular  data  set.  The  drawing  of  the  axes  and  zero  reference 
lines  is  handled  in  the  same  manner  as  in  DRAW1  (see  sect.  2.1.1). 

The  addition  ot  a border  to  the  graph  and  the  use  of 
different  colored  ink  are  handled  by  the  same  formal  parameter.  The 
two  options  can  be  selected  individually  or  together.  When  axes  are 
drawn  in  the  standard  position,  a border  is  drawn  at  the  top  and  right 
of  the  graph  connected  to  the  major  axes.  This  border  consists  of  axes 
of  the  same  type  (linear  or  log)  as  those  directly  opposite,  but 
without  annotation.  When  the  option  to  change  pen  color  is  selected,  a 
PM  message  is  printed  on  the  operator's  console  after  the  graph  is 
labeled  and  the  axes  drawn,  and  before  each  data  set  is  plotted.  The 
PM  message  puts  the  output  stream  in  a wait  status,  the  text  of  the 
message  requests  the  operator  to  change  the  pen  on  the  plotter  (the 
user  must  have  told  the  operator  in  advance  which  colors  to  use  in  what 
order) . There  is  no  automatic  pen  selector  with  this  plotter.  When 
the  graph  is  completed,  another  PM  message  requests  the  restoration  of 
the  original  color. 

Finally,  scaling  factors  are  controlled  through  the  third 
section.  These  are  used  to  normalize  the  data  set(s)  being  plotted 
(this  scaling  factor  is  the  power  of  10  that  appears  as  part  of  the 
axis  label) , to  set  the  increments  that  determine  the  scale  for  a 
linear  axis,  and  to  set  the  minimum  and  maximum  powers  of  10  for  a 
logarithmic  axis.  These  four  values  must  be  encoded  into  one  integer 
word  ("INCR"  in  this  software)  of  twelve  digits.  Each  value  is 
allotted  three  digits.  For  linear  axes,  the  factor  for  the  independent 
axis  occupies  the  high-order  (left-most)  digits,  followed  by  the 
increment  for  the  independent  scale,  the  scaling  factor  for  the 
dependent  axis,  and  the  increment  for  the  dependent  scale.  For 
logarithmic  axes,  the  minimum  power  of  10  for  the  independent  axis 
occupies  the  high-order  digits,  followed  by  the  maximum  power  of  10 
for  the  independent  axis,  the  minimum  power  of  10  for  the  dependent 
axis,  and  the  maximum  power  of  10  for  the  dependent  axis.  The  range 
that  the  exponents  for  the  scaling  factors  for  the  log  scale  can  cover 
is  from  +99  to  -99.  To  cover  the  positive  range,  the  high-order  digit 
of  the  appropriate  set  of  three  digits  must  be  zero.  For  the  negative 
range,  the  high-order  digit  of  the  appropriate  set  of  three  digits  must 
be  one.  An  exponent  of  zero  is  achieved  by  setting  the  appropriate 
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digits  to  200  (setting  the  digits  to  zero  causes  the  INCR  to  be 
ignored) . For  a logarithmic  scale  the  uses  of  a partial  decade 
preceding  the  minimum  value  or  following  the  maximum  value  are  not 
considered  if  999  is  encoded  into  the  appropriate  portion  of  INCR.  If 
the  data  set(s)  being  plotted  cannot  meet  the  user-defined  scaling 
f actor (s),  scale  increment (s) , or  minimum  and  maximum  scale  values,  the 
supplied  values  are  ignored  and  new  values  are  computed.  A printed 
message  notifies  the  user  of  this  action.  The  values  used  to  define 
the  scaling  factors  are  as  follows: 


exponent  of  10  (10  ) to  which  the  independent 

arrays  are  normalized 

scale  increment  for  the  independent  axis 
POWY 

exponent  of  10  (10  ) to  which  the  dependent  arrays 

are  normalized 

scale  increment  for  the  dependent  axis 
minimum  power  of  10  for  the  independent  axis 

maximum  power  of  10  for  the  independent  axis 

minimum  power  of  10  for  the  dependent  axis 

maximum  power  of  10  for  the  dependent  axis 


Order  of  encoding  of  information: 


POWX  AX  POWY  AY  linear- linear  plot 
MINX  MAXX  MI NY  MAXY  log-log  plot 
MINX  MAXX  POWY  AY  log-linear  plot 


POWX  AX  MINY  MAXY  linear-log  plot 


Specific  examples  of  the  defining  of  scaling  factors  follow. 

(a)  To  have  the  independent  array (s)  normalized  by  10"6  and  a 
scale  increment  of  20,  with  the  dependent  array (s)  normalized  by  103  and 
a scale  increment  of  5, 


POWX  will  equal 

106  (100  indicates  negative  power 

10 

of 

AX  will  equal 

020 

POWY  will  equal 

003 

AY  will  equal 

005 

INCR  will  equal 

106020003005 

(b)  To  have  the  independent  array  (s)  normalized  by  10°  with  no 
specified  scale  increment,  and  the  dependent  array (s)  plotted  with  DRAW 4 
software  calculated  values. 

POWX  will  equal 

200  (200  indicates  a factor  of  1) 

AX  will  equal 

000 

POWY  will  equal 

000 

AY  will  equal 

000 

• 

• • 

INCR  will  equal 

200000000000 

(c)  To  have  independent  array (s)  plotted  on  a logarithmic  scale 
covering  from  10" 1 to  104,  and  dependent  array Cs)  plotted  on  a 

logarithmic  scale  covering  from  10'6  to  10’1  , 

MINX  will  equal 

101 

MAXX  will  equal 

004 

MINY  will  equal 

106 

MAXY  will  equal 

101 

INCR  will  equal 

101004106101 

15 


(d)  To  have  the  independent  array (s)  plotted  on  a linear  scale 
normalized  by  10* 6 with  no  specified  scale  increment  and  the  dependent 
array (s)  plotted  on  a logarithmic  scale  starting  at  101  with  no 
specified  maximum  and  no  partial  decade  with  the  calculated  maximum. 


POWX 

will 

equal 

106 

will 

equal 

000 

MINY 

will 

equal 

001 

MAXY 

will 

equal 

999 

INCR 

will 

equal 

106000001999 

The  values  contained  in  the  scratch  arrays  provided  to 
DRAW4  in  the  third  section  (IP=3)  are  restored  if  the  value  for  the 
number  of  points  in  the  scratch  arrays  is  positive.  This  is  done  by 
storing  the  values  in  the  scratch  arrays  on  the  same  scratch  disc  file 
used  to  store  previous  plot  information,  and  restoring  the  values  to  the 
scratch  arrays  before  returning  to  the  calling  program  from  DRAW4.  If 
the  value  for  the  number  of  points  in  the  scratch  arrays  is  negative, 
the  values  for  the  scratch  arrays  are  not  restored. 

The  FORTRAN  subprogram  DRAW 4 is  called  as  follows: 

For  IP=1 : 


CALL  DRAW4 ( IP , IFILE , IX LAB , IYLAB , IPTLAB , ISBLAB , XLAB , YLAB , 
PTLAB, SUBLAB) 


where 


IFILE  A number  from  1 to  8 identifying  the 

scratch  file  associated  with  the  graph. 
For  each  graph  the  value  must  be 
consistent  for  each  value  of  IP. 

IXLAB  Number  of  words  in  the  label  for  the 

independent  axis,  not  to  exceed  a value  of 
3.  If  IXLAB=0 , axis  label  is  omitted. 

IYLAB  Number  of  words  in  the  label  for  the 

dependent  axis,  not  to  exceed  a value  of 
3.  If  IYLAB=0,  axis  label  is  omitted. 
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IPTLAB 


number  of  words  in  the  plot  title,  not  to 
exceed  a value  of  8.  If  IPTLAB=0,  plot 
title  is  omitted. 


ISBLAB  Number  of  words  in  the  plot  subtitle,  not 
to  exceed  a value  of  8.  If  ISBLAB=0,  plot 
subtitle  is  omitted. 

XLAB  Array  containing  user-defined  label  for 

the  independent  axis,  not  to  exceed  3 
words  (30  characters) . If  IXLAB=0,  a 
dummy  name  can  be  used  in  the  call. 

YLAB  Array  containing  user-defined  label  for 

the  dependent  axis,  not  to  exceed  3 words 
(30  characters).  If  IYLAB=0,  a dummy  name 
can  be  used  in  the  call. 

PTLAB  Array  containing  user-defined  title  for 

the  plot,  not  to  exceed  8 words  (80 
characters) . If  IPTLAB=0,  a dummy  name 
can  be  used  in  the  call. 

SUBLAB  Array  containing  user-defined  subtitle  for 
the  plot,  not  to  exceed  8 words  (80 
characters) . If  ISBLAB=0,  a dummy  name 
can  be  used  in  the  call. 

For  IP=2 : 

CALL  DRAW4 (IP,IFILE,ITYPE, NPTS , KODE , I SPACE , XVAL , 

YVAL , YMIN , YMAX ) 

where 

IFILE  A number  from  1 to  8 identifying  the 

scratch  file  associated  with  the  graph  and 
a previous  call  to  DRAW4  with  IP=1. 

I TYPE  Code  used  to  select  the  type  of  plot. 

Read  only  on  the  first  call  with  IP=2  for 
each  IFILE. 

If  ITYPE=1  a linear  plot  is  selected. 

=2  a log-log  plot  is  selected. 
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=3  a linear-log  plot  (semilog) 
with  the  linear  axis  for  the 
dependent  array)  is  selected. 

=4  a log-linear  plot  (semilog 
with  the  linear  axis  for  the 
independent  array)  is 

selected. 

NPTS  Number  of  points  in  both  the  independent  and 
dependent  data  arrays  to  be  plotted.  For 
multiple  plots , NPTS  for  each  data  set  does 
not  have  to  be  the  same. 

KODE  The  character  to  be  drawn  at  data  points, 
connected  by  a line.  If  KODE  is  negative, 
only  the  sumbol  is  drawn. 


If  KODE=0 

no  character 

=1 

a plus  (+) 

=2 

an  X 

=3 

a triangle  (A) 

=4 

a square  (□) 

=5 

an  hour  glass  (X) 

=6 

an  up  arrow  (t) 

=7 

a star  (*) 

Other  characters  can  be  drawn  by  equating  KODE  to  a 
character  (e.g. , 2RJJA,  where  'A'  is  the  character 
to  be  drawn  at  data  points) . 


ISPACE  (a)  ISPACE/10  is  the  rate  of  occurrence  for  the 
drawing  of  the  symbol  designated  by  the 
variable  KODE.  If  ISPACE/10>0  and  KODE=0, 
intermediate  points  are  disregarded  to  give 
the  user  a simplified  plot. 
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(b)  The  last  digit  of  ISPACE  determines  the 
type  of  line  to  be  drawn,  provided  KODE  is 
not  negative. 

If  the  last  digit  of  ISPACE  is 

0 a solid  line  is  drawn  

1 a dashed  line  is  drawn  

2 a dotted  line  is  drawn.  . . . 

3 a dash-dot  line  is  drawn  

XV AL  Array  containing  values  of  independent  variables 

to  be  plotted. 

YVAL  Array  containing  values  of  dependent  variables 

to  be  plotted. 

YMIN  Used  in  the  first  call  with  IP=2  to  set  the 

minimum  value  for  the  dependent  axis.  If  YMIN 
is  zero  (0. ) it  is  ignored.  For  a minimum  value 
of  zero,  set  YMIN  equal  to  a very  small  number 
(e.g. , l.E-30) . Any  values  in  the  arrays  that 
are  less  than  YMIN  are  set  equal  to  YMIN,  thus 
modifying  the  arrays. 

YMAX  Used  in  the  first  call  with  IP=2  to  set  the 

maximum  value  for  the  dependent  axis.  If  YMAX 
is  zero  (0.)  it  is  ignored.  Any  values  in  the 
arrays  that  are  greater  than  YMAX  are  set  equal 
to  YMAX,  thus  modifying  the  arrays. 

For  IP=3 : 

CALL  DRAW4 (IP , IFILE, IDUM, INCR, ICHOICE ,N,X,Y ,ZAX ,ZAY) 
where 


IFILE  A number  from  1 to  8 identifying  the  scratch 
file  associated  with  the  graph  and  previous 
calls  to  DRAW 4 with  IP=1  and  2. 

IDUM  Dummy  variable,  presently  not  used. 


19 


\ 


INCR 


ICHOICE 


N 

X 

Y 

ZAX 


Variable  used  to  specify  scales  other  than  those 
calculated  by  the  subprogram.  A single  integer 
number  composed  of  four  values  (POWX,  AX,  POWY, 
AY,  MINX,  MAXX,  MINY,  or  MAXY)  which  are  used 
for  the  normalizing  factors  and  scale 
increments;  ignored  if  equal  to  zero  (0). 

Variable  used  to  request  a change  in  pen  color 
and  to  draw  borders  on  graph. 


If  ICHOICE=0  no  action  is  taken. 


=1 

a change  of 

pen 

is  requested 

before  the  plotting 

of  each  data 

set . 

=2 

plot  is  drawn 

with  a border 

(upper  and 

possible. 

right 

sides) , if 

=3 

both  options 
implemented. 

1 

and  2 are 

Number  of  points  in  scratch  arrays;  should  be 
equal  to  or  greater  than  largest  data  set  being 
plotted.  If  N is  negative  the  scratch  arrays 
are  not  restored  to  their  original  values. 

Scratch  array  used  for  temporary  storage  of 
independent  array (s)  to  be  plotted. 

Scratch  array  used  for  temporary  storage  of 
dependent  array (s)  to  be  plotted. 

Variable  that  determines  where  the  independent 
axis  is  drawn. 

If  ZAX=0.  axis  is  drawn  in  the  standard 
location.  This  is  also  the  default 
position  when  the  next  two  conditions 
cannot  be  met. 

=1.  axis  is  drawn  at  the  zero  occurrence 
of  the  dependent  axis. 
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* 


ZAY 


=2.  axis  is  drawn  at  standard  location 
with  the  addition  of  a zero  reference 
line  placed  at  the  zero  occurrence  on 
the  dependent  axis. 

Variable  that  determines  where  the  dependent 
axis  is  drawn. 

If  ZAY=0.  axis  is  drawn  in  standard  location. 

This  is  also  the  default  position 
when  the  next  conditions  cannot  be 
met . 


=1 . axis  is  drawn  at  the  zero 
occurrence  of  the  independent  axis. 

=2.  axis  is  drawn  at  standard  location 
with  the  addition  of  a zero 
reference  line  placed  at  the  zero 
occurrence  on  the  independent  axis. 


The  following  messages  cam  be  printed  to  the  line  printer  by  DRAW4 : 

(a)  "CALLS  OUT  OF  SEQUENCE  ON  FILE  XX" 

This  message  is  printed  if  the  second  section  of  DRAW4  is  called 
before  the  first  section,  or  the  third  section  called  before  the  first 
or  second  sections,  for  the  graph  associated  with  scratch  file  number 
XX. 

(b)  "XX  POINTS  DROPPED  FROM  TRACE  YY  IN  FILE  ZZ , NEG.  OR  ZERO 

VALUE  IN  ARRAY" 

This  message  notifies  the  user  of  XX  invalid  values  (negative  or 
zero)  for  the  plotting  of  the  YYth  data  set  on  a log  scale  for  the 
graph  associated  with  scratch  file  number  ZZ. 

(c)  "TRACE  ON  FILE  XX  AFTER  TRACE  NO.  YY  CANNOT  BE  PLOTTED, 

NEG.  OR  ZERO  VALUE  IN  ARRAY" 

This  message  is  printed  when  an  invalid  value  (negative  or  zero) 
occurs  after  valid  values  in  data  set  YY+1  on  a log  scale  for  the  graph 
associated  with  scratch  file  number  XX. 

(d)  "AMIN=XX  AMAX=YY" 

This  message  provides  the  user  with  additional  information  on  the 
minimum  and  maximum  values  of  an  array  when  their  difference  is  very 
small  and  expansion  of  the  resolution  of  the  scale  would  cause  the 
scale  notation  to  overlap  illegibly. 
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(e)  "IMPROPER  SCALE  IGNORED’ 


This  message  notifies  the  user  that  the  selected 
scale  could  not  be  used.  The  program  then  attempts 
scale  close  to  the  sele9ted  values. 


values  for  the 
to  calculate  a 


(f)  "ARRAY  PASSED  IN  CALL  3,  FILE  XX,  NOT  LARGE  ENOUGH" 

This  message  warns  the  user  that  the  dimension  of  the  scratch 
arrays  provided  in  the  third  section  associated  with  file  number  XX  is 
not  as  large  as  that  of  the  data  sets  being  plotted.  The  results  are 
unpredictable . 

2.1.3  Subroutine  DRAW2 

Subroutine  DRAW23  is  a printer-plot  version  of  DRAW1.  The 
printer-plot  routine  mimics  DRAW1  and  is  useful  for  making  test  runs  or 
plots  for  which  lower  resolution  is  acceptable.  The  calling  parameters 
for  the  printer-plot  routine  are  identical  to  those  for  DRAW1  and  the 
routine  is  accessed  by  replacing  the  calls  to  DRAW1  with  calls  to 
DRAW  2. 

2.1.4  Subroutine  DRAW 3 

Subroutine  DRAW33  is  a printer-plot  version  of  DRAW4.  The 
printer-plot  routine  mimics  DRAW4  and  is  useful  for  making  test  runs  or 
plots  for  which  a lower  resolution  is  acceptable.  The  calling 
parameters  for  the  printer-plot  routine  are  identical  to  those  for 
DRAW 4 and  the  routine  is  accessed  by  replacing  the  calls  to  DRAW 4 with 
calls  to  DRAW 3. 

2.1.5  Subroutine  CAL CMP 

Subroutine  CALCMP  translates  calls  to  "CALCMP"  (the 
California  Computer  Products,  Inc.,  plot  subroutine)  into  a form 
compatible  with  the  Houston  Instrument  plotter.  Image  scaling  and  beam 
control  are  not  possible,  and  reference  to  these  functions  causes  a 
stop  condition.  To  implement  this  routine  the  user  simply  loads  from 
the  ANAPAC  library  before  loading  from  the  CALCOMP  library. 


3Egon  Marx,  Printer  Version  of  Plots  Made  by  an  Incremental  Plotter, 
Harry  Diamond  Laboratories  TM-75-33  (December  1975). 


22 


2.2  Fourier  Transform  Routines 


2.2.1  Subroutine  FFT 


Subroutine  FFT  performs  a forward  or  an  inverse  Fourier 
transform  of  an  equispaced  array  of  N points,  where  N=2**NP0W. 
Subroutine  FFT  makes  use  of  the  Cooley-Tukey  algorithm,  programmed  in 
COMPASS  assembly  language.  The  input  and  output  for  FFT  is  passed  by 
the  same  formal  parameter,  A.  For  the  inverse  transform,  the  output  is 
obtained  by  taking  the  REAL  part  of  the  output  array.  By  convention, 
the  sign  of  the  exponent  in  the  integral  is  negative  for  the  forward 
transform  and  positive  for  the  inverse  transform. 

Subroutine  FFT  is  called  as  follows: 

CALL  FFT (A, NPOW.N, DX , ISIGN) 

A Complex  array  containing  input  or  output  data. 

NPOW  Power  of  2. 

N Number  of  points  in  data  array. 

DX  Increment  for  independent  array. 

ISIGN  Plus  or  minus  1;  value  of  p desired  in  transform, 
/ g (t)  eP^i-f  fcdt . 

2.2.2  Subroutine  FLAT 


Subroutine  FLAT4  performs  a forward  Fourier  transform  of  an 
equispaced  array  by  approximating  the  given  function  with  a piecewise 
linear  function  and  then  using  the  Cooley-Tukey  algorithm  to  obtain  a 
sampling  of  the  Fourier  transform.  The  input  and  output  for  FLAT  are 
passed  by  the  same  formal  parameter,  A. 

Subroutine  FLAT  is  called  as  follows: 

CALL  FLAT ( A , N , DX , IS IGN ) 

A Complex  array  containing  input  and  output 

data. 

N Number  of  points  in  data  array;  a power  of 

2. 

4 Alfred  Brandstein  and  Egon  Marx,  Numerical  Fourier  Transform,  Harry 
Diamond  Laboratories  TR-1748  (December  1975). 
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DX  Increment  for  independent  array. 

ISIGN  Plus  or  minus  1;  value  of  p desired  in 
transform  /g  (t)  eP^if t<3t . 

2.2.3  Subroutine  FLIT 

Subroutine  FLIT4  performs  an  inverse  Fourier  transform  of  an 
equispaced  array.  The  input  and  output  for  FLIT  are  passed  by  the  same 
formal  parameter,  A.  The  input  data  are  placed  in  the  first  N/2+1 

points.  The  output  is  obtained  by  taking  the  REAL  part  of  the  whole 
output  array. 

Subroutine  FLIT  is  called  as  follows: 


CALL  FLIT(A,N,DX, ISIGN) 


A 

Complex  array  containing  data. 

N' 

Number  of 

points  in  data  array;  a 

power  of  2. 

DX 

Increment 

for  independent  array. 

ISIGN 

Plus  or 

minus  1;  value  of  p 

desired  in 

transform 

, /g(f)eP2iriftdt. 

2.2.4  Subroutines  NUFT  and  INUFT 

Subroutines  NUFT  and  INUFT 4 perform  a direct  and  an  inverse 
Fourier  transform  of  a function  defined  by  straight  lines  connecting 
the  data  points.  Equispaced  arrays  and  a number  of  points  which  is  a 
power  of  2 are  not  a prerequisite  of  either  subroutine. 

Subroutine  NUFT  is  called  as  follows : 

CALL  NUFT (X, Y , N1 ,N2 , OM, IOM , FT, DX, JFLAG) 

X Independent  array  containing  time  data, 

input . 

Y Dependent  array  containing  amplitude  data, 

input . 


4 Alfred  Brandstein  and  Egon  Marx,  Numerical  Fourier  Transform,  Harry 
Diamond  Laboratories  TR-1748  (December  1975) . 
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N1  Number  of  points  in  input  arrays. 

N2  Number  of  point  in  output  arrays. 

OM  Independent  array  containing  frequency  data, 
input . 

IOM  Variable  that  indicates  whether  OM  array  is 
frequency  or  circular  frequency  for  input  and 
output. 

If  I0M=1  circular  frequency  input  and 

output . 

=2  circular  frequency  input, 

frequency  output. 

=3  frequency  input,  circular 

frequency  output. 

=4  frequency  input  and  output. 

FT  Complex  array  for  transformed  data,  output. 

DX  Scratch  array,  dimensioned  to  Nl. 

JFLAG  Plus  or  minus  1,  sign  of  exponent  for 
transform. 

Subroutine  INUFT  is  called  as  follows: 

CALL  INUFT (X,Y,Nl,N2,0M, IOM , FT , DX , JFLAG) 

X Independent  array  for  time  data,  output. 

Y Dependent  array  for  amplitude  data,  output. 

Nl  Number  of  points  in  output  arrays. 

N2  Number  of  points  in  input  arrays. 

OM  Independent  array  containing  frequency  data, 
input . 

IOM  Variable  that  indicates  whether  OM  array  is 
frequency  or  circular  frequency  for  input  and 
output. 
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frequency  input  and 


FT 


If  I0M=1 

circular 

output. 

=2 

circular 

frequency 

=3 

frequency 

frequency 

=4 

frequency 

frequency  input, 
output. 

input,  circular 

ouptut . 

input  and  output. 


Complex  array  containing  transformed  data, 
input. 


DX  Complex  scratch  array,  dimensioned  to  N2. 

JFLAG  Plus  or  minus  1,  sign  of  exponent  of 
transformed  data. 


2.2.5  Subroutine  FFTA 


Subroutine  FFTA  performs  a forward  Fourier  transform  or  an 
inverse  Fourier  transform  with  full  or  half  aliasing.  For  a forward 
transform,  FFTA  simply  calls  FFT  with  a negative  exponent.  For  half 
aliasing,  the  values  in  the  second  half  of  the  transform  array,  on  the 
midpoint  plus  one,  are  replaced  with  the  conjugate  of  the  respective 
value  of  the  first  half.  For  full  aliasing,  the  conjugates  of  the 
values  in  the  first  half  of  the  transform  array,  on  the  midpoint  plus 
one,  are  added  to  the  respective  values  in  the  second  half  of  the 
transform  array.  Then  the  conjugates  of  the  sums  for  each  value  in 
the  second  half  replace  the  respective  values  in  the  first  half.  Full 
aliasing  should  be  used  only  on  analytic  functions.  The  result  of 
either  the  half  or  the  full  aliasing  is  then  used  as  input  to  FFT.  A 
phenomenon  that  can  occur  with  inverse  FFT's  is  an  overall  amplitude 
shift  of  the  resulting  function.  To  adjust  for  this  shift,  the  formal 
parameter  ZRO  allows  the  user  to  shift  the  total  function  so  that  the 
zero  time  point  has  an  amplitude  of  zero. 

Subroutine  FFTA  is  called  as  follows : 

CALL  FFTA(K,FT,N,DX, ISIGN, ZRO,Y) 

K Variable  that  selects  the  type  of  transform 

and  type  of  aliasing. 
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If  K=0 


forward  Fourier  transform  selected. 

=1  inverse  Fourier  transform  with  full 
aliasing  selected. 

=2  inverse  Fourier  transform  with  half 
aliasing  selected. 

FT  Complex  array  of  time  data  input  and 

transform  data  output. 

N Number  of  points  in  input  and  output  arrays. 

DX  Increment  for  independent  arrays. 

ISIGN  Plus  or  minus  1,  sign  of  exponent  of 

transformed  data. 

ZRO  Logical  variable  used  to  select  option  to 

shift  output  time  function. 

If  ZRO=.T.  function  shifted  so  that  first 
point  has  amplitude  of  zero. 

=.F.  function  not  shifted. 

Y Output  array  containing  real  part  of  inverse 

Four;-  transform. 

2.2.6  Subroutine  SMU7 

The  power  spectrxxa  from  a digitized  time-amplitude  trace 
shows  a strong  oscillating  noise  component,  especially  on  a logarithmic 
scale.  Subroutine  SMUZ  can  be  used  to  present  the  output  in  a more 

intelligible  form.  The  input  is  an  array  of  a function  given  at 

constant  intervals.  When  two  maxima  or  cwo  minima  occur  closer  than  a 
prescribed  number  of  points  (the  threshold) , the  function  in  between  is 
replaced  by  an  average  value  between  the  straight  line  joining  the 
maxima  and  that  joining  the  minima,  with  a tapered  beginning  and  end. 
The  process  is  repeated  until  there  are  no  changes  in  a complete  pass. 
Two  difference  thresholds  can  be  prescribed  for  two  sectors  of  the 
function.  The  number  of  passes  through  the  procedure  is  printed  in  the 
output  (usually  four  to  six  passes) . The  thresholds  have  to  be  chosen 
in  such  a way  that  the  unwanted  noise  is  eliminated  while  the 
significant  extrema  remain;  a number  to  start  with  might  be  1/100  of 

the  total  number  of  points.  It  has  been  found  that  the  use  of  this 

subroutine  on  the  real  and  imaginary  parts  of  a Fourier  transform  that 
subsequently  has  to  be  inverted  tends  to  introduce  spurious 
oscillations  for  late  times.4 

4 Alfred  Brandstein  and  Egon  Marx,  Numerical  Fourier  Transform,  Harry 
Diamond  Laboratories  TR-1748  (December  1975). 
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Subroutine  SMUZ  is  called  as  follows : 


r 


; 


: 


: 


CALL  SMUZ(X,N,M1,M2,FRAC) 

X Ordinates  of  function  being  smoothed. 

N Number  of  points  in  X array. 

Ml  Threshold  number  of  points  for  first  part  of 

curve. 

M2  Threshold  number  of  points  for  second  part 

of  curve. 

FRAC  Fraction  of  points  in  first  part  of  curve. 
2.2.7  Subroutine  OMEGA 

Subroutine  OMEGA  generates  a frequency  array.  Given  the 
maximum  frequency  desired,  WMAX,  and  the  initial  delta  frequency,  DW, 
OMEGA  constructs  an  unequispaced  array  with  increments  in  geometric 
progression.  If  DW  equals  zero,  an  equispaced  frequency  array  is 
constructed.  The  parameters  WMAX  and  DW  can  be  in  either  units  of 
frequency  or  circular  frequency. 

Subroutine  OMEGA  is  called  as  follows: 

CALL  OMEGA (OM, WMAX, N,DW) 

OM  Output  frequency  array. 

WMAX  Maximum  frequency  desired. 

N Number  of  frequency  points  desired. 

DW  Initial  frequency  step. 

2.3  Data  Preparation  and  Interpolation  Routines 
2.3.1  Subroutine  CSTOUT 

Subroutine  CSTOUT  checks  for  time  ordering  of  the  independent 
array  and  casts  out  those  points  not  in  an  ascending  time  order.  If 
points  are  encountered  with  the  same  time  value,  the  corresponding 
values  in  the  dependent  array  are  averaged.  Messages  notify  the  user 
how  many  points  were  valid  and  invalid.  Subroutine  CSTOUT  is 
recommended  for  use  with  digitized  data. 
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Subroutine  CSTOUT  is  called  as  follows: 


CALL  CSTOUT (X,Y,N) 

X Independent  array,  time. 

Y Dependent  array,  amplitude. 

N Number  of  points  in  X and  Y arrays;  if  invalid 
values  encountered,  new  value  of  N returned. 

2.3.2  Subroutines  CLINT  and  LINT 

Subroutines  CLINT  and  LINT  perform  a linear  interpolation 
between  the  points  in  the  given  array.  The  outputs  from  CLINT  and  LINT 
are  equispaced  arrays.  The  output  array  for  CLINT  is  complex;  that  for 
LINT  is  real.  Any  points  which  fall  outside  the  range  of  the  input 
variable  are  given  a value  of  zero  in  the  output  array. 

Subroutines  CLINT  and  LINT  are  called  as  follows: 

CALL  CLINT (X,Y ,N,M,DX,YOUT) 


or 

CALL  LINT(X,Y ,N,M,DX,YOUT) 

X Independent  array,  time. 

Y Dependent  array,  amplitude. 

N Number  of  points  in  X and  Y arrays. 

M Number  of  points  in  YOUT  array. 

DX  Increment  for  independent  variable. 

YOUT  Dependent  output  array;  complex  for  CLINT. 

2.4  File  Structuring  Routines 

The  following  subroutines  provide  a compatible  data  structure 
for  the  easy  exchange  of  data  between  different  computer  programs. 
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Subroutine  WRTAPE  writes  data  sets  onto  a local  program  file 
which  is  to  be  established  as  a permanent  file  for  access  by  another 
program.  Multiple  data  sets  can  be  written  on  the  same  file.  The  user 
must  properly  close  the  file  before  termination  of  the  program. 

Subroutine  WRTAPE  is  called  as  follows: 

CALL  WRTAPE ( NT, X,Y,N, LABEL) 

NT  File  number  where  data  are  written. 

X Independent  array;  for  an  equispaced  array 

X(l)  contains  the  increment. 

Y Dependent  array. 

N Number  of  points  in  data  arrays;  negative 

value  for  N indicates  an  equispaced 
independent  array. 

LABEL  Name  of  data  set  (eight  words) . 

2.4.2  Subroutine  RDTAPE 

Subroutine  RDTAPE  reads  data  sets  from  the  permanent  files 
created  by  using  WRTAPE.  Subroutine  RDTAPE  checks  for  the  end  of 
information  (EOI)  and  for  any  irregularities  in  the  data  format.  This 
subroutine  reads  digitized  data  files. 

Subroutine  RDTAPE  is  called  as  follows: 

CALL  RDTAPE (NT, X,Y,N, LABEL) 

NT  File  number  where  data  are  contained. 

X Independent  array;  for  equispaced  array  X(l) 

contains  the  increment. 

Y Dependent  array . 

N Number  of  points  in  data  arrays;  if  N is 

negative,  independent  array  is  equispaced. 

LABEL  Name  of  data  set  (eight  words) . 
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2.5 


Analytic  Signal  Generating  Routines 


A model  used  to  represent  analytically  the  response  to  EMP  is 
a function  that  vanishes  for  t < 0,  defined  by 


f(t)  = At/Tj 


0 < T < Tj 


'n(t  - Tj ) I 

= Aexp  t — C l (t  - TjJJ  cos  2(T2  - Ti~)J  Tl£  5 T2 


= a + 6t  + <5t2  + yt3 


T2<t<  T3 


= aexp[-C2(t  - T3)]  cosfwjtt  - T3)] 


+ (a£2/w2)  exp[-53(t  - T3)]  sin[u2(t  - T3)] 


T3  < t < te 


=0 


t > T 


The  values  of  A,  Tlf  T2,  T3,  and  a are  defined  by  the  user.  The 
user  also  gives  a time  T to  define  the  end  of  the  "recorded"  part  of 
the  pulse.  The  values  of  E£lf  C2  and  £3  are  given  through  constants  Sj, 
S2  and  S3  by 

r _ 1 ■ - 

1 " Si(T2  - T2) 


52  s2(TE-T3) 


s3(TE  - t3)  ' 
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and  the  circular  frequencies  are  given  by 


Wl  ~ 2<T \ - T3) 


rrS4 


Finally,  c»,  6,  y,  and  6 are  determined  by  matching  the  value  and  the 
derivative  of  the  function  at  the  ends  of  the  interval. 

The  Fourier  transform  of  this  function  is  the  sum  of  the  Fourier 
transforms  of  the  components  in  the  different  intervals.  These 
components  are 


Fj  (w) 


A[exp(-iu)Tj)  (iwTj  +1)  - 1] 


TjU)2 


F2  (U>) 


Aexp(-iu)Tj) 


< 


(Ci  + iw) ' 


x {exp[- (Ci  + iw)  (T2  - Ti)]  2CT2^~-~r'i)  + (Sl  + iu)  * 


(u) 


exp(-iuiT3) 


2C  + 6yT3 


io) 


3 


- exp(iioT2) 


26  + 6yT2 
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F 4(10)  = aexp(-id)T3) 


[C2  + iw 
(C2  + iu)2  + o»i‘ 


e2  exp[-(£2  + iw)  (T  - T3I2 

+r)  V It — 


(£3  + id))2  + w22  (£2  + ioi)z  + o>i 


4 2 x 


X [-(c2  + iw)  cosdij  (Te  - T3)  + ui  sind)j  (te  - Tj|] 

£2  exp^-(C3  + id))  (Te  - T3)j 
d>2  (£3  + id))2  + d)22 

x[-(£3  + id))  sind)2(TE  - T3j  - d>2  coso)2  (T£  - T3)]j  , 


where 


ttA 


2 (T2  - Tj)  exPl-^l(T2  - T1>1  ' 


2a  Y 

Y + 


(t3  - T2)  2 (t3  - T2)  3 


6 = 


(T3  - T2)“ 


- Y(T2  + 2T3)  , 


and  n is  1 or  0 depending  on  whether  we  assume  that  the  trace  vanishes 

after  T or  is  given  by  the  expression  in  the  last  interval. 

K 
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r 


1 


Functions  Fi  and  F3  have  to  be  computed  separately  for  w = 0. 
They  become 


ATi 

Fl<0)  =— 


6 (T32  - T22)  (t32  - t23)  6(t34-t24) 

F 3 (0)  =a  (t3  - T2)  + + - + 


where 


6 = -2T36  - 3T32y  , 


a = -T2B  - t225  - T23y. 


The  function  is  continuous  at  the  extremes  of  the  intervals, 
but  it  does  not  necessarily  vanish  at  T . The  derivative  is 
discontinuous  at  Tj,  but  continuous  at  T2  (where  the  function  vanishes) 
and  at  T3  (where  the  function  reaches  a minimum  given  by  a) . At  T4,  the 
first  term  in  the  corresponding  expression  for  f(t)  vanishes. 

2.5.1  Subroutine  ANMI 


Subroutine  ANMI4  computes  the  value  of  f(t)  for  the 
analytical  model  described  in  section  2.5.  The  desired  time  points  and 
the  number  of  points  and  input  to  ANMI  and  the  amplitude  values  are 
returned.  Subroutine  ANMI  requires  a NAMELIST  card,  PULSE,  that 
contains  the  values  A,  AA,  Tl,  T2,  T3,  T4,  TE,  SI,  S2,  S3,  and  S4.  The 
meanings  of  these  are  evident,  except  that  AA  = a. 

Subroutine  ANMI  is  called  as  follows: 

CALL  ANMI (X, Y,N) 

X Time  points  at  which  function  evaluated. 

Y Array  for  amplitudes. 

N Number  of  points  in  X. 


4 Alfred  Brandstein  and  Egon  Marx,  Numerical  Fourier  Transform,  Harry 
Diamond  Laboratories  TR-1748  (December  1975) . 


2.5.2  Subroutine  ANTRA 


Subroutine  ANTRA4  computes  the  Fourier  transform  for  the 
analytical  model  described  in  section  2.5.  The  desired  frequency  points 
(circular  frequency)  and  the  number  of  points  are  input  to  ANTRA  and  the 
values  of  the  Fourier  transform  are  returned.  Subroutine  ANTRA  requires 
a NAMELIST  card,  PULSE,  which  contains  the  values  A,  AA,  Tl,  T2,  T3,  T4, 
TE,  SI,  S2,  S3,  S4,  and  FINIT.  Variable  FINIT  takes  the  value  0.  for  an 
infinite  trace  and  1.  for  a finite  trace. 

Subroutine  ANTRA  is  called  as  follows: 

CALL  ANTRA (OM, FT, N) 

OM  Circular  frequencies  at  which  function 
evaluated. 

FT  Complex  array  containing  values  of  the 
Fourier  transform. 

N Number  of  points  in  OM. 

3.  THE  TRANL  USER'S  LIBRARY  AND  SUPPORTING  PROGRAMS 

The  TRANL  user's  library  is  a collection  of  independent  computer 
programs  with  supporting  subroutines  which  provides  the  means  for  the 
validation  of  digitized  data  sets  and  the  means  to  easily  incorporate 
digitized  data  into  the  normal  program  input  structure.  The  library 
comprises  three  principal  programs.  The  main  program,  TRANS,  translates 
the  punched  card  output  from  the  digitizing  system,  checks  for  format 
errors  in  the  data  set,  provides  a time-amplitude  listing  of  the 
digitized  data,  and  creates  a local  file  of  the  digitized  data  which  is 
formatted  to  comply  with  subroutines  WRTAPE  and  RDTAPE  in  the  ANAPAC 
library.  If  errors  are  encountered  during  the  translation  of  a 
digitized  data  set,  the  user  is  informed  of  the  card  number  where  the 
error  occurred.  To  aid  in  the  correction  of  the  error  the  program  GRAFC 
can  be  used  to  translate  the  digitized  data  set  into  the  control  format 
used  in  the  digitization  process.  Another  supporting  program,  TEDD, 
provides  for  "time  window  expansion  of  digitized  data."  These  programs 
were  stored  in  TRANL  with  an  access  level  of  one,  which  allows  access  to 
the  programs  by  control  cards. 

Another  program  associated  with  the  TRANL  library,  but  not  stored 
in  it,  is  REDUCE.  Program  REDUCE  provides  the  user  with  the  ability  to 
selectively  reduce,  by  defining  time  windows,  the  number  of  points  in  a 
digitized  data  set  and  create  a new  data  file. 


4 Alfred  Brandstein  and  Egon  Marx,  Numerical  Fourier  Transform,  Harry 
Diamond  Laboratories  TR-1748  (December  197 5). 


35 


V 


3.1  Program  TRANS 

Program  TRANS  translates  the  punched  card  output  from  the 
digitizing  system  into  time-amplitude  information.  The  time -amplitude 
information  is  listed  on  the  printer  and  written  on  a local  file,  which 
can  be  catalogued  as  a permanent  file.  The  printed  listing  can  be 
eliminated  by  putting  sense  switch  1 to  "on."  The  local  or  permanent 
file  can  then  be  accessed  by  the  use  of  RDTAPE  from  ANAPAC.  The  input 
for  TRANS  can  be  one  or  more  data  sets  each  preceded  by  a label  card  and 
followed  by  an  end-of-information  card  (EOI,  7/8/9  punch).  The  input 
stream  is  terminated  by  two  EOI  cards. 

Program  TRANS  is  accessed  by  control  cards  as  follows: 

TRANS (lfnl , lfn2 , lfn3) 

lfnl  Logical  file  name  for  input  (e.g.,  INPUT). 

Ifn2  Logical  file  name  for  output  (e.g.,  OUTPUT). 

Ifn3  Logical  file  name  assigned  by  user  to  receive 
time-amplitude  data. 

3.2  Program  GRAFC 

Program  GRAFC  converts  the  punched-card  digitized  data  into 
the  control  commands  and  data  points  recorded  by  the  digitizer.  The 
output  listing  from  GRAFC  greatly  aids  in  error  corrections  and 
modification  of  digitized  data  and  can  eliminate  the  need  for 
re-digitization.  The  input  to  GRAFC  is  the  same  as  for  TRANS. 

Program  GRAFC  is  accessed  by  control  cards  as  follows: 

GRAFC  (lfnl,lfn2) 

lfnl  Logical  file  name  for  input  (e.g.,  INPUT). 

Ifn2  Logical  file  name  for  output  (e.g.,  OUTPUT). 

3.3  Program  TEDD 

Program  TEDD  provides  the  user  with  the  ability  to  plot 
selected  time  windows  of  a digitized  data  set.  The  user  can  thus  expand 
portions  of  a time-amplitude  trace  for  clarity  or  display  the  individual 
digitized  points  in  the  selected  time  window.  By  displaying  individual 
digitized  points,  invalid  points  can  be  identified  and  either  be  removed 
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from  the  data  deck  or  have  their  values  changed.  To  aid  in  the 
correction  of  a digitized  trace,  TEDD  prints  the  following  information 
about  a data  set: 


(a)  The  label  assigned  to  the  data  set 

(b)  The  number  of  digitized  points  in  the  trace 

(c)  The  duration  of  the  trace 


(d)  The  selected  boundary  points  for  the  time  window 

(e)  The  time  of  the  end  points  in  the  digitized  trace  which 
bound  the  selected  time  window  and  define  the  time  scale 
for  the  plotted  output 

(f)  The  number  of  digitized  points  in  the  time  window 

tg)  The  numbered  location  of  the  digitized  end  points  relative 
to  the  beginning  of  the  trace 

Program  TEDD  can  process  multiple  data  sets  with  a maximum  of 
six  time  windows  per  data  set  (the  time  windows  can  overlap) . The 
maximum  and  minimum  for  the  dependent  variable  can  also  be  specified  to 
maintain  the  same  perspective  for  the  plots.  The  format  for  the  maximum 
and  minimum  values  and  the  time  windows  is  ' 6X,2E10.3'  with  one  time 
window  per  card.  Each  group  of  time  windows  must  be  separated  by  an  EOI 
card  (7/8/9  punch).  By  putting  sense  switch  2 to  "on,"  the  marking  of 
the  individual  data  points  is  eliminated. 

The  input  data  for  TEDD  are  as  follows : 

Card  1 AMIN  Minimum  value  of  scale  for  plotting  of 
dependent  array;  ignored  if  equal  to  zero. 

AMAX  Maximum  value  of  scale  for  plotting  of 
dependent  array;  ignored  if  equal  to  zero. 

Format  (6X,2E10.3) 

Card  2 T1  Lower  limit  of  first  time  window  to  be 

expanded . 


T2  Upper  limit  of  first  time  window  to  be 

expanded . 

Format  (6X,2E10.3) 
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Card  3 


T1&T2 


time 


window; 


same  as 


Limits  for  second 
Card  2. 


Card  7 T1&T2  Limits  for  sixth  time  window;  same  as 

Card  2. 

Program  TEDD  can  also  be  applied  to  other  types  of  data  files 
compatible  with  RDTAPE  from  the  ANAPAC  library.  Program  TEDD  requires 
the  support  of  the  ANAPAC  library. 

Program  TEDD  is  accessed  by  control  cards  as  follows: 

TEDD (lfnl , If n2 , lfn3) 

lfnl  Logical  file  name  for  input  (e.g.,  INPUT). 

Ifn2  Logical  file  name  for  output  (e.g.,  OUTPUT). 

Ifn3  Logical  file  name  assigned  by  user  which  contains 
data  to  be  examined. 

3.4  Program  REDUCE 

Program  REDUCE  enables  the  user  to  selectively  reduce  the 
number  of  points  in  a digitized  record.  Program  REDUCE  operates  in  two 
modes.  The  first  mode  allows  the  user  to  define  from  one  to  four  time 

windows  and  the  percentage  of  data  points  to  be  saved  from  each  time 

window.  In  the  second  mode,  the  user  provides  the  factor  by  which  the 

trace  is  to  be  reduced.  The  reduced  time  trace  is  then  written  on  a 

user-defined  logical  file  using  WRTAPE  and  also  plotted  to  display  the 
results  of  the  reduction.  Program  REDUCE  requires  the  support  of  the 
ANAPAC  library. 

The  input  data  for  REDUCE  are  as  follows: 

Card  1 Factor  Factor  by  which  data 

set  is  to  be  reduced;  if 
0 , time  windows  and 
percentages  defined  on 
following  cards. 

M Number  of  points  to  which 

data  set  is  to  be  reduced; 
can  be  omitted  if 
Factor^O . 

Format  (I2,2X,I4) 
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Card  2 (if  Factor=0)  Ti  Upper  limit  of  time 


window;  i=l,2,3;  last  time 
window  bounded  by  Ti  and 
maximum  time  of  data  set 
if  EPi^lOO  percent. 

Format  (3E10.3) 


Card  3 (if  Factor=0)  Pi  Percentage  of  points 

selected  from  corre- 
sponding time  windows  of 
i;  i=l,2,3;  for  each  Ti 
there  must  be  a Pi. 

Format  (3E10.3) 


Abnormal  termination  occurs  under  the  following  conditions; 


STOP  1 First  percentage  value  equal  to  zero. 


STOP  2 Ti  exceeds  the  r.^xim\om  time  value  of  the  data 
set. 


Program  REDUCE  is  stored  as  a permarent  file  and  is  accessed 
by  a Scope  ATTACH  card.  Loading  and  execution,  are  accomplished  with  a 
Scope  file  name  call.  Access  is  accomplished  as  follows: 

ATTACH (A , REDUCE , id=***** , MR=1 ) 

A ( , , lfnl , lfn2) 


lfnl  Logical  file  name  assigned  by  user  which  contains 
data  set  to  be  reduced. 


Ifn2  Logical  file  name  assigned  by  user  which  receives 
reduced  data  set. 


4.  USE  OF  THE  ANAPAC  AND  TRANL  LIBRARIES  BY  GAMPHS 


1 


The  computer  analysis  section  of  the  GAMPHS1  application  (see 
fig.  3)  makes  use  of  the  TRANL  user's  library  and  many  of  the 
subroutines  from  the  ANAPAC  user's  library.  The  routines  in  the  TRANL 
user's  library  (see  sect.  3)  and  the  digitizing  system  that  the  routines 
support  provide  the  link  between  system  field  tests  and  the  computer 
analysis  of  the  field  test  data.  The  output  of  the  TRANL  routines  is  a 
file  of  digital  data.  The  ANAPAC  user  library  (see  sect.  2)  provides 
many  of  the  analytic  routines  required  to  analyze  the  test  data.  The 
modular  ANAPAC  routines  are  used  as  building  blocks  to  develop  programs 
which  yield  meaningful  information  about  the  data  (e.g.,  Fourier 
transform,  transfer  functions) . The  routines  also  prepare  the  data  to 
serve  as  input  for  further  analysis  by  other  computer  programs.  The 
file  structuring  routines  contained  in  the  ANAPAC  user's  library  (see 
sect.  2.4)  provide  the  file  compatibility  necessary  to  pass  data  between 
computer  programs. 


Figure  3.  Interaction  of  computer  programs  used  by  GAMPHS. 


1 George  Gornak  et  al,  EMP  Assessment  for  Army  Tactical  Communi- 
cations Systems:  Transmission  Systems  Series  No.  1 — Radio  Terminal  Set 

AN/TRC-145  (U) , Harry  Diamond  Laboratories  TR-1746  ( February  1976). 
(SECRET-RESTRICTED  DATA) 
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The  following  are  other  computer  programs  which  are  part  of  GAMPHS 
and  use  portions  of  the  TRANL  and  ANAPAC  user's  libraries  for 
analytic  support  and  file  compatability : 

1.  GROUND:1 * * * 5  Determines  the  total  electric  and  magnetic  fields 

in  air  and  those  transmitted  below  the  surface  for 
a plane-wave  pulse  incident  on  a ground  of  finite 
conductivity.  The  output  can  be  used  as  input  for 
TEMPO  or  FREFLD.  Uses  Fourier  transform,  plotting 
and  file  structuring  routines. 

2.  TEMPO:1  Computes  the  transient  response  of  a cable, 

antenna,  or  aperture  to  EMP  stimulation;  both  a 
batch  and  interactive  code.  The  output  can  be  used 
as  input  for  DAMTRAC . Uses  plotting  and  file 
structuring  routines. 

3.  FREFLD:6  Calculates  the  current  response  of  shielded  and 

unshielded  cables  to  EMP  stimulation.  The  output 
can  be  used  as  input  for  DAMTRAC.  Uses  Fourier 
transform,  plotting  and  file  structuring  routines. 

4 . DAMTRAC : 7 Performs  damage  analysis  of  modeled  solid-state 

circuits  from  EMP  stimulation.  Uses 
data -preparation,  plotting  and  file-structuring 
routines.  See  appendix  A of  this  report  for  the 
use  of  digitized  data  by  DAMTRAC. 


1 George  Gornak  et  al,  EMP  Assessment  for  Army  Tactical  Communi- 

cations Systems : Transmission  Systems  Series  No.  1 — Radio  Terminal  Set 

AN/TRC-145  (U) , Harry  Diamond  Laboratories  TR-1746  (February  1976). 

(SECRET-RESTRICTED  DATA) 

5 Egon  Marx,  Reflected  and  Transmitted  Fields  for  a Plane-Wave  Pulse 
Incident  on  a Conducting  Ground,  Harry  Diamond  Laboratories  TR-1740 
(April  1975) . 

6 Robert  F.  Gray,  Nuclear  Electromagnetic  Pulse  Simulation  by  Point 
Source  Injection  Technique  for  Shielded  and  Unshielded  Penetration, 
Harry  Diamond  Laboratories  TR-1737  (December  1975). 

7 George  Baker  et  al , Damage  Analysis  Modified  TRAC,  Harry  Diamond 
Laboratories  TM-76-6  (May  1975). 
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APPENDIX  A. “COMPATIBILITY  OF  DAMTRAC  WITH  DIGITIZED  DATA  FILES 
AND  THE  PROPER  APPLICATION  OF  THE  DATA 


To  enhance  user  ability  to  describe  transient  circuit  driving 
functions,  DAMTRAC1  has  been  designed  to  accept  time-domain  data  files 
structured  for  WRTAPE/RDTAPE  compatibility  (see  sect.  2.4  of  the  body  of 
the  report) . Time-domain  data  files  compatible  with  the  WRTAPE/RDTAPE 
format  include  digitized  data  files  and  data  files  generated  by  the 
programs  covered  by  the  GAMPHS  application  (see  sect.  4 of  the  body  of 
the  report) . During  a circuit  analysis  run,  the  digitized  or  foreign 
data  file  is  "attached"  and  copied  to  logical  unit  7.  The  DAMTRAC 
subroutine  DGSRCE  reads  the  time-amplitude  data  from  logical  unit  7 and 
edits  the  data,  reducing  the  number  of  points  to  500,  if  the  input  data 
exceed  500.  However,  the  editing  scheme  used  by  DGSRCE  is  to  delete 
every  other  point  for  files  containing  between  500  and  1000  points  and 
to  delete  every  third  point  for  files  containing  more  than  1000  points. 
This  method  of  reducing  the  number  of  points  used  for  the  analysis  can 
inadvertently  alter  the  signature  of  the  transient  driving  function. 
For  this  reason  the  use  of  the  program  REDUCE  (see  sect.  3.4  of  the  body 
of  the  report) , which  selectively  reduces  the  number  of  time-amplitude 
points  and  displays  the  results,  is  recommended  for  the  reduction  of 
driving  function  inputs  to  DAMTRAC. 


1 George  Baker  et  al , Damage  Analysis  Modified  TRAC,  Harry  Diamond 
Laboratories  TM-76-6  (May  1975). 
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ARLINGTON,  VA  20360 
ATTN  PME-117-21,  SANGUINE  DIV 

HEADQUARTERS,  NAVAL  MATERIEL  COMMAND 
STRATEGIC  SYSTEMS  PROJECTS  OFFICE 
1931  JEFFERSON  DAVIS  HIGHWAY 
ARLINGTON,  VA  20390 
ATTN  NSP2201,  LAUNCHING  & HANDLING 

BRANCH,  BR  ENGINEER,  P.  R.  FAUROT 
ATTN  NSP-230,  FIRE  CONTROL  & GUIDANCE 
BRANCH,  BR  ENGINEER,  D.  GOLD 
ATTN  NSP-2701,  MISSILE  BRANCH, 

BR  ENGINEER,  J.  W.  PITSENBERGER 


COMMANDER 

HQ  SPACE  AND  MISSILE  SYSTEMS  ORGANIZATION 
P.O.  96960  WORLDWAYS  POSTAL  CENTER 
LOS  ANGELES,  CA  90009 
ATTN  S7H , DEFENSE  SYSTEMS  APL  SPO 
ATTN  XRT,  STRATEGIC  SYSTEMS  DIV 
ATTN  SYS,  SURVIVABILITY  OFC 

SPACE  AND  MISSILE  SYSTEMS  ORGANIZATION 
NORTON  AFB,  CA  92409 
ATTN  MMH,  HARD  ROCK  SILO  DEVELOPMENT 

COMMANDER 

AF  SPECIAL  WEAPONS  CENTER,  AFSC 
KIRTLAND  AFB,  NM  87117 


HARRY 

ATTN 


ATTN 

ATTN 

ATTN 

ATTN 

ATTN 

ATTN 

ATTN 

ATTN 

ATTN 

ATTN 

ATTN 

ATTN 


DIAMOND  LABORATORIES 
MCGREGOR,  THOMAS,  COL,  COMMANDING 
OFFICER/FLYER,  I. N. /LANDIS,  P.E./ 
SOMMER,  H. /CONRAD,  E.E. 

CARTER,  W.W.,  DR.,  ACTING  TECHNICAL 
DIRECTOR/MARCUS,  S.M. 

KIMMEL,  S.  , IO 
CHIEF,  0021 
CHIEF,  0022 
CHIEF,  LAB  100 
CHIEF,  LAB  200 
CHIEF,  LAB  300 
CHIEF,  LAB  400 
CHIEF,  LAB  500 
CHIEF,  LAB  600 
CHIEF,  DIV  700 
CHIEF,  DIV  800 
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DISTRIBUTION  (Cont'd) 


HARRY  DIAMOND  LABORATORIES  (Cont'd) 

ATTN  CHIEF,  LAB  900 
ATTN  RECORD  COPY,  BR  041 
ATTN  HDL  LIBRARY  (3  COPIES) 

ATTN  CHAIRMAN,  EDITORIAL  COMMITTEE 

ATTN  CHIEF,  047 

ATTN  TECH  REPORTS,  013 

ATTN  PATENT  LAW  BRANCH,  071 

ATTN  MCLAUGHLIN,  P.W. , 741 

ATTN  LANHAM,  C.,  PROGRAM  & PLANS  OFFICE 

ATTN  CHIEF,  0024 

ATTN  CHIEF,  1020  (20  COPIES) 

ATTN  CHIEF,  1030 
ATTN  CHIEF,  1040 
ATTN  CHIEF,  1050 
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